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Aims. We study the physical characteristics (shape, dimensions, spin axis direction, albedo maps, mineralogy) of the dwarf-planet 
Ceres based on high-angular resolution near-infrared observations. 

Methods. We analyze adaptive optics J/H/K imaging observations of Ceres performed at Keck II Observatory in September 2002 with 
an equivalent spatial resolution of ~50 km. The spectral behavior of the main geological features present on Ceres is compared with 
laboratory samples. 

Results. Ceres' shape can be described by an oblate spheroid (a = b = 479.7 ± 2.3 km, c = 444.4 ±2.1 km) with EQJ2000.0 spin 



vector coordinates a = 288° ± 5° and So = +66° ± 5°. Ceres sidereal period is measured to be 9.074 10+jj 



h. We image surface 



features with diameters in the 50-180 km range and an albedo contrast of ~6% with respect to the average Ceres albedo. The spectral 
behavior of the brightest regions on Ceres is consistent with phyllosilicates and carbonate compounds. Darker isolated regions could 
be related to the presence of frost. 
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1. Introduction 

Ceres is by far the largest body among the population of main- 
belt asteroids. Curiously, although it was discovered more than 
200 years ago [Piazzi, 1802], significant progress to understand 
its surface and interior properties has only been made over the 
last two decades. In particular we have seen a renewal of in- 
terest towards Ceres triggered by (i) the availability of sensi- 
tive spectro-imaging instrumentation on medium-to-large size 
telescope, enabling a detailed study of Ceres surface and physi- 
cal properties [e.g. Mitchell et al., 1996; Drummond et al., 1998; 
Dotto et al., 2000, and others], and (ii) more recently the selec- 
tion of the NASA DAWN Discovery mission, which will visit 
Vesta and Ceres in 2011 and 2015 respectively [Russell et al., 
2003]. But the main scientific interest drawn by Ceres is that it 
provides an excellent laboratory to understand how planetoids 
accreted early in the history of our solar system, and the role 
of volatiles in planetary formation and evolution. Unlike Vesta, 
which is dry and shows evidence of the melting phase and plan- 
etary differentiation [Keil, 2002] seen in larger terrestrial plan- 
ets like Earth, the "dwarf planet" Ceres (as it should now be 
named following the IAU guidelines) shows strong signs of wa- 
ter alteration on its surface [Jones et al., 1990]. In this respect, 
Ceres displays stronger similarities with the icy outer satellites 
of Jupiter than with the dry asteroids that populate the inner re- 
gion of the Main Belt. A possible scenario is that Ceres formed 
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in a "wet" environment, from the accretion of both rocky plan- 
etesimals originally present at this heliocentric distance in the 
early planetary nebula, and icy planetesimals that migrated in- 
ward from the outer regions and whose ices had been preserved 
[Mousis & Alibert, 2005]. As a consequence of this "wet" his- 
tory, the study of Ceres is of paramount importance to under- 
stand the process of planetary accretion and formation of the 
low-albedo primitive asteroids that populate the outer part of the 
Main Belt. 

With a typical angular diameter of 0.6" at opposition, Ceres 
can be spatially resolved from the ground using adaptive optics 
instruments available on medium-to-large telescopes. We thus 
carried out a program of multi-band imagery of Ceres from Keck 
Observatory at high-angular resolution, with the goals of pre- 
cisely deriving its shape, dimensions, direction of spin axis, and 
distributions of albedo and color across its surface. 

2. Observations and data reduction 

2.1. Observations 

Our Ceres observations were made with the Keck II telescope 
on Mauna Kea, Hawaii, on 2002 September 22 and 28 UT, 
one week before the opposition of Ceres. These observations 
were obtained under optimal atmospheric conditions with a 
sub-arcsec seeing at an airmass lower than 1.65 (with half of 
the data taken with an airmass lower than 1.2). Whereas the 
2002 opposition occured near Ceres aphelia (Ceres was at a 
geocentric and heliocentric distance of 1.98 AU and 2.94 AU 
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Filter 


A c 


AA 





Across 


0Over 




G"m) 


(pm) 


(km) 


Diameter 


Surface 


J 


1.248 


0.163 


37.2 


-26 


-666 


H 


1.633 


0.296 


47.9 


-20 


-400 


K 


2.124 


0.351 


64.6 


-15 


-228 



Table 1. Central wavelength (A c ) and bandpass width (AA) for 
each filter. The equivalent size (in km) of the theoretical reso- 
lution element (0) on Ceres and the number of resolution el- 
ements across the diameter and over the apparent disk of Ceres 
are also reported. Ceres covers more than 3 200 pixels (projected 
major- and minor-radius of -33 and -31 pixels respectively) on 
the NIRC2 detector. 



respectively), its apparent angular diameter was -0.66". The 
phase angle was of 7° and 5.5° for the two nights, leading 
to an illuminated fraction of the surface of 99.6% and 99.7% 
respectively. 

Ceres was imaged using NIRC2, the second-generation 
near-infrared camera (1024x1024 InSb Aladdin-3) and the 
adaptive optics (AO) system installed at the Nasmyth focus of 
the Keck II telescope [van Dam et al., 2004]. The images of 
Ceres were acquired at 3 near-infrared wavebands J [1.166- 
1.330 fan], H [1.485-1.781 fan], and K [1.948-2.299 fan], with 
an image scale of 9.942 + 0.050 milliarcsec per pixel. Within 
the Ceres observation sequence we interspersed observations of 
reference stars, at similar airmass and through the same set of 
filters, to evaluate the system Point Spread Function (PSF). This 
calibration was required to perform a posteriori image restora- 
tion as described in the next section. No offset to sky was done 
but the telescope position was dithered between each exposure 
in order to record simultaneous sky and object frames, while the 
target (science or calibration) was positioned at three different 
locations on the detector, separated by nearly 5" from each other. 



2.2. Data reduction 

We first reduced the data using the standard procedure for near- 
infrared images. A bad pixel mask was made by combining the 
hot and dead pixels found from the dark and flatfield frames. The 
bad pixels in our calibration and science images were then cor- 
rected by replacing their values with the median of the neighbor- 
ing pixels (7x7 pixel box). Our sky frames were obtained from 
the median of each series of dithered science image, and then 
subtracted from the corresponding science images to remove the 
sky and instrumental background. By doing so, the dark current 
was also removed. Finally, each image was divided by a normal- 
ized flatfield to correct the pixel-to-pixel sensitivity differences 
of the detector. 

After these first basic reduction steps, we applied image 
deconvolution techniques to our set of Ceres data using the 
Mistral algorithm [Fusco, 2000; Mugnier et al., 2004]. The use 
of such an algorithm permits to restore the optimal spatial reso- 
lution of each image and is particularly well adapted to decon- 
volve objects with sharp edges, such as asteroids. Image restora- 
tion techniques are known to be constrained by the limitations to 
measure the precise instrumental + atmosphere responses at the 
exact same time the science observations are made. Mistral is 
a myopic deconvolution method, which estimates both the most 
probable object, and the PSF, from analysis of science and refer- 
ence star images. The time needed for the algorithm to converge 
is largely dependent on the image size. Due to our large num- 




Fig. 1. Selected views of Ceres in K-band during -200° of its 
rotation. The image is oriented with North up and East left. 
The values of the Sub-Earth Point longitude (SEP^) measured 
positively from to 360 degrees in a right-hand system (fol- 
lowing the IAU recommandation [Seidelmann et al., 2005]) are 
indicated below their corresponding images. Our meridian ori- 
gin is chosen to be the same than Li et al. [2006a]. The Ceres 
spin axis and positive pole are also indicated. The image stretch 
was chosen to enhance the surface features visible on Ceres at 
the detriment of the terminator. The two main surface features 
present on Ceres are indicated with the arrows A and B (see 3.1) 
and can be followed during part of their rotation. The brighter 
spots visible near the limb of Ceres in some of the images are 
artifacts from the deconvolution (see 4.2). 



ber of images to process, we decreased the deconvolution com- 
putation time by resizing all our images to a smaller (128x128 
pixels, but still over 2x the diameter of Ceres) window centered 
on the object (Ceres or PSF). Comparison tests showed that no 
deconvolution artifacts were introduced by the use of smaller 
images. The deconvolved images of Ceres were then compared 
to identify the few frames (9 in total, that is less than 2.5% of 
the full set of data) whose outputs were not satisfactory and dis- 
card them from our set of data. In the end we obtained 360 im- 
ages of Ceres with a spatial resolution (Table 1) equivalent to 
the diffraction limit of a 10 m telescope (given by the angular 
sampling © = A/D, with A the wavelength and D the telescope 
diameter). A subset of the restored images is presented in Fig. 1. 
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Fig. 2. Apparent motion of the two main features A (bright) and 
B (dark) over ~90° of Ceres rotation during the night of the 2002 
September 22 UT (SEP^ a +8°, p„ 338°). We used the rota- 
tional track of these two features to determine the direction of 
the rotation axis (as described in text). 



3. Ceres global physical properties 

3.1. Spin vector coordinates 

Measuring the Spin Vector Coordinates (SVC) of Ceres has 
always been a difficult task. The small amplitude of its 
lightcurve (~0.04 mag) prevented its determination using classi- 
cal lightcurve inversion techniques. Until recently, the solutions 
that were reported for the pole coordinates of Ceres were widely 
dispersed (spanning a 90° range in the plane of the sky during 
the 1995 opposition as reported by Parker et al. [2002], based 
on the compilation of Johnson et al. [1983], Millis et al. [1987], 
Saint-Pe et al. [1993] and Drummond et al. [1998] pole solu- 
tions). This situation changed dramatically with the use of high- 
contrast direct imaging techniques capable of resolving the disk 
of Ceres and image fine details on its surface, such as AO obser- 
vations (e.g [Dumas et al., 2003], and this paper) and HST ob- 
servations (e.g. [Parker et al., 2002; Thomas et al., 2005]). The 
high-quality data obtained at Keck allow us to obtain an inde- 
pendent measurement of the SVC of Ceres, which is compared 
to the HST results [Thomas et al., 2005]. Precise knowledge of 
the SVC is mandatory to re-project the individual images into 
Ceres planetocentric referential and establish albedo maps of its 
surface (see 4.2). 

We performed a control-point method, which consists of fol- 
lowing the apparent motion of a surface feature while the as- 
teroid rotates. The path of a rotating feature on the surface of 
Ceres describes an ellipse if the asteroid is seen under a Sub- 
Earth Point latitude (SEP^) different from 0°, or ±90° (where it 
describes a line, or a circle respectively). We marked the posi- 
tions of two main features A and B (Fig. 1) while they rotated 
respectively over a 93.9° and 87.7° range around Ceres' spin 
axis (100 and 87 images respectively). All these measurements 
are shown together on a 3-D representation of Ceres in Fig. 2. 
The shape of the projected track being determined by the direc- 
tion of the rotational axis, we were able to measure the north 
pole angle (p„) and the SEP^ at the time of the observations. 

In order to compare our solution with other results, we fol- 
lowed the IAU recommendations and calculated the SVC in the 
J2000.0 equatorial frame: a = 288° + 5° and 6 = +66° ± 5°. 
This result is in agreement with the latest result from HST 
[Thomas et al., 2005] (a = 291° ± 5° and 6 = +59° + 5°). 



Epoch 


Date 
(Julian Day) 


SEP, 

O 




ti 


2452 539.89402 


234.6 ± 


2 


t 2 


2452 545.915 46 


261.4 ± 


2 


t 3 


2453 002.24128 


22.6 ± 


5 



Table 2. Sub-Earth Point longitude (SEP )) for the three epochs 
(ti and t2 from Keck, t3 from HST). Using ti as reference, Ceres 
sidereal phasing with ?2 and ?3 can be obtained by applying a 
+ 1.3° and a -98.7° correction taking into account the apparent 
geometry of the Earth with respect to Ceres. 



Sidereal period 


SEP,(t 2 ) 


ASEP, 


SEP,(t 3 ) 


ASEP, 


(h) 


C) 


O 


o 


o 


9.066 588 


255.1 


6.3 


17.9 


4.7 


9.066 685 


255.1 


6.3 


22.6 


0.0 


9.066 780 


255.1 


6.3 


27.2 


-4.6 


9.074000 


259.8 


1.6 


17.8 


4.8 


9.074090 


259.8 


1.6 


22.1 


0.5 


9.074 100 


259.8 


1.6 


22.6 


0.0 


9.074110 


259.8 


1.6 


23.1 


0.5 


9.074200 


259.9 


1.6 


27.5 


-4.9 


9.081526 


264.4 


-3.0 


17.9 


4.7 


9.081526 


264.5 


-3.1 


22.6 


0.0 


9.081526 


264.6 


-3.2 


27.2 


-4.6 



Table 3. Sidereal periods of Ceres which are in agreement with 
our observations, based on SEP^fi) = 234.6°. The predicted 
Sub-Earth Point longitude (SEP^) are computed at epochs t2 and 
t3 and are compared with the measurements (difference ASEP.i). 
We only kept sidereal periods whose predicted SEP i(t) were in- 
side measurement uncertainties. 



To directly derive the obliquity of a planetary body, it is also 
convenient to express the SVC in an ecliptic reference frame. 
Our measurements give Aq = 7° and /3o = +83° (with 5° un- 
certainty), leading to an obliquity for Ceres of about 4°. Such 
a small obliquity was expected from analysis of lightcurve data 
[Tedesco et al., 1983] and thermal properties of Ceres regolith 
[Spencer, 1990]. Indeed, a larger obliquity would imply stronger 
variations in the amplitude of Ceres' lightcurves, which have 
never been reported. As a result, Ceres can only be observed over 
a small range of Sub-Earth Point latitude (SEP^) [-11°, +10°] 
(computed until 2015), imposing on us to wait for the arrival of 
the DAWN spacecraft to get a direct view of its polar regions. 

3.2. Rotational period 

We establish a precise and independent measurement of the ro- 
tation period of Ceres using images from our two Keck epochs 
(2002 September 22 and 28 UT), plus one additional processed 
Hubble ACS/HRC image (j8p502amqJofMt) from the HST pro- 
gram GO 9748 retrieved from the Small Body Nodes archive 
[Li et al., 2006b] for the epoch 2003 December 28 UT (Table 2) 
[Thomas et al., 2005; Li et al., 2006a; Parker et al., 2006]. Two 
main albedo marks visible in all images were used to obtain pre- 
cise information on Ceres rotational phase at these three epochs 
(see Table 2). In addition, the period between the Keck and HST 
observations was sufficiently large to accurately measure the er- 
ror accumulation on Ceres' period over more than 1000 rota- 
tions. 

We then used the Eproc ephemeris generator [Berthier, 
1998] to predict the longitude of the Sub-Earth Point and de- 
termine the value of the Ceres sidereal period that minimizes the 
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a 


6 


W 


Ps 


Epoch 


C) 


C) 


C) 


(h) 


(JD) 


288 ±5 


+66 ±5 


-46 ±2 


9.0741 ± 10~ 4 


2452 539.894 02 



Table 4. Ceres sidereal period (P s ), pole direction (ao, So) and 
initial rotational phase angle (Wo) at epoch of reference in the 
J2000.0 equatorial frame. 



difference ASEP^ between the observed and the computed SEP,}. 
We adopted the pole solution derived in this work as well as the 
priorly determined Ceres sidereal period of P s = 9.075 + 10~ 3 
h [Tedesco et al., 1983] as a best estimate. Table 3 shows the 
sidereal periods that are in agreement with the Sub-Earth Point 
longitude (SEP^) at the epochs of the observations. 

The value of Ceres sidereal period best matching our obser- 
vations was found to be P s = 9.074 10+°-°°°}° h with a ~ 10~ 4 
h resolution. The small ~0.5 s uncertainty, which come mainly 
from the error on the mesured longitudes (see Table 2), makes 
possible to predict Ceres' Sub-Earth Point longitude, SEP,}, with 
an error of only 40° in rotational phase over the next decade. By 
comparison, the very recent study by Chamberlain et al. [2007] 
compiled 50 years of lightcurve measurements to derive a period 
of P s = 9.074 170+2 x 10^ 6 h. 



3.3. Dimensions 

Precise measurements of the shape of Ceres allow us to remotely 
investigate its internal structure and test whether it is a differenti- 
ated body [Thomas et al., 2005]. We thus performed a Laplacian 
of Gaussian (LoG) wavelet analysis to 192 images of Ceres to 
extract its limb contours. Because it is very sensitive to varia- 
tions of gradient, the LoG permits to precisely detect the inflec- 
tion points in the flux distribution of our individual deconvolved 
images of Ceres. We found Ceres to be rotationally symmetric as 
first reported from a preliminary analysis of the Keck- AO data 
by Dumas et al. [2003], and later confirmed by Thomas et al. 
[2005]. Analysis of our complete set of contours did not return 
any deviation from our ellipsoidal model larger than ~18 km 
(see Fig. 3). From this model, and our knowledge of the direc- 
tion of its spin axis, we were able to determine the minor- and 
major- projected radius of Ceres and correct them from their as- 
pect and phase angles. We found that Ceres is well described by 
an oblate spheroid whose semi-axes are a = b — 479.7 ± 2.3 km 
andc — 444.4+2.1 km (l-<x dispersion for the fitted axis). These 
values are different from HST's by ~ 1 km (relative difference of 
two percent). In comparison, the agreement between our semi- 
major axis and the determination made from stellar occultation 
by Millis et al. [1987] is remarkable: 479.6 ± 2.4 km (occulta- 
tion) vs 479.7 + 2.3 km (Keck), while our minor-axis value differ 
from theirs by -6-9 km (Table 5). Whereas Millis et al. [1987] 
assumed a zero-obliquity at the time of the occultation, our SVC 
solution gives a SEP^, of +3.3°. The corrected polar radius from 
stellar occultation is thus 453.3 + 4.5 km, which is still differ- 
ent from our measurement by nearly 10 km. It is important to 
note that during the occultation of 1987, the cords along Ceres 
orbit were aligned in an East- West direction, whereas the north 
pole position angle was p n ~ 339°. This implies that the estima- 
tion of the small radius of Ceres was more loosely constrained 
than its semi-major axis. In addition, Table 5 provides a com- 
parison between the pixel size of the HST and Keck images, and 
shows that the coarser sampling of the disk of Ceres as seen by 
HST could explain the differences in our respective determina- 




50 100 150 200 250 300 
Sub-Earth Longitude 



Fig. 3. Left figure: Extracted limb contour (black) from a J-band 
image of Ceres with its ellipsoidal fit overplotted for compari- 
son (dotted gray). This example shows that deviations from the 
ellipsoidal model are much smaller than the resolution element 
(gray disk inside the box) obtained with Keck-AO. Right fig- 
ure: Compilation of our deviation measurements as function of 
the Sub-Earth Point longitude (SEP,}). No deviations larger than 
~18 km (half a resolution element at J-band) are detected. 



Quantity 


Ima 


gery 


Occultation 


(km) 


Keck 


HST 


Sol.l 


Sol.2 


a 


479.7 ± 2.3 


487.3 ± 1.8 


479.6 ± 2.4 


481.6 ±2.4 


c 


444.4 ±2.1 


454.7 ± 1.6 


453.4 ± 4.5 


450.1 ± 2.0 


R 


467.6 ± 2.2 


476.2 ± 1.7 


470.7 ± 3.1 


470.8 ± 2.3 


pixel 


14.3 ±0.7 


33.9 x 29.6 







Table 5. Semi-major (a) and -minor (c) axes, and equiv- 
alent radius (R- yaac) for Ceres derived in this paper 
(Keck), Thomas et al. [2005] (HST), and Millis et al. [1987] 
(Occultation). The pixel size of the Keck and HST images are 
also given for comparison. 



tions of its size. Our mean radius for Ceres is estimated to be 
R = 467.6 ± 2.2 km, which is also closer to the determination 
made from stellar occultation than HST's. 

Adopting a mass for Ceres ofM = 9.43 + 0.07 x 10 20 kg (av- 
erage of most recent measurements [Viateau & Rapaport, 1998; 
Michalak, 2000; Kovacevic & Kuzmanoski, 2007]), we find a 
mean density p = 2 206 + 43 kg.irT 3 . This value is relatively 
high for a hydrated G-type asteroid like Ceres, but can be ex- 
plained by a low porosity [see Britt et al., 2002], and is simi- 
lar to the density of the icy outer Jovian satellites Ganymede 
and Callisto. One can assume Ceres to be in hydrostatic equi- 
librium and inverse the relation between a, c and J 2 [given by 
McCord & Sotin, 2005] as following 



h = 



1 - 



27T 2 fl 3 



a PlGM 



2a 



and find J 2 = 26.7 x 10 3 + 1.9 x 10 5 . By comparison, 
Thomas et al. [2005] have found J 2 = 21.7x 10^ + 8.5 x 10~ 5 . If 
we refer to McCord & Sotin [2005], these two independent es- 
timates of J 2 correspond to their internal models #2 and #3 for 
Ceres, that is to a differentiated Ceres with a silicate-rich region 
in its center. The rather large difference between the two J 2 de- 
terminations is due to the fact that J 2 is highly dependant on the 
a/c ratio which shows large uncertainties depending on the tech- 
nique used for its measurement. Lightcurve analysis for a quasi- 
spherical object is generally little sensitive to this parameter; 
while measurements made from direct imaging from the Earth 
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distance suffer from lack of sufficient spatial resolution to esti- 
mate this parameter with the required precision. Nevertheless, 
both studies converge towards a differentiated asteroid and grav- 
ity field measurement made by DAWN will provide a better un- 
derstanding of its mass repartition and internal structure. 

Conrad et al. [2007] have shown that a detailed study of an 
asteroid shape, and its departure from a perfect ellipsoid, can be 
carried out using images deconvolved with Mistral. Deviations 
of Ceres limb measurements from our shape model can thus di- 
rectly be linked to topography. Fig. 3 shows that no features de- 
viating from our shape model by more than 15-18 km can been 
observed. The highest relief expected on Ceres is calculated to 
be ~ 10-20 km high [see also Johnson & McGetchin, 1973, Fig. 
1], thus confirming that no significant deviation from an ellip- 
soid can be detected given the size of our resolution element. 
Search for relief would require a resolution of at least about 5 
km to provide unambiguous detection of a topography. 

4. Ceres surface composition 

4.1. Background 

Ceres occupies a particular place in our solar system. It is 
physically located far enough from the Sun to have been 
preserved from strong heating during the T Tauri phase of 
the Sun [Ghosh et al., 2006, and references therein], and has 
possibly retained some of its primordial elements. Its low am- 
plitude lightcurve (~0.04 mag), which cannot result from Ceres 
symmetrical shape [Thomas et al., 2005], shows that its surface, 
while rather uniform, displays faint albedo features, unlike the 
igneous asteroid Vesta whose hemispheric albedo variations 
are among the strongest seen among main belt objects [Gaffey, 
1997]. Similarly, in contrast to Vesta, whose shape is highly 
irregular due to past collisions [Thomas et al., 1997], Ceres 
displays a uniform spheroidal shape, deprived of strong surface 
and topographic features, and its density (see 3.3) cannot simply 
be attributed to macroporosity of its internal material [see 
Britt et al., 2002]. All these characteristics point to the presence 
of volatile elements in the interior of Ceres. Mousis & Alibert 
[2005] show that Ceres could have accreted from an assemblage 
of rocky and icy planetesimals, even at such short heliocentric 
distances. They calculated that icy planetesimals could have 
drifted from more distant regions of our solar system (up to 
15 AU heliocentric) to the actual position of Ceres (2.7 AU) 
without losing entirely their volatiles. This idea supports the 
model of Fanale & Salvail [1989] in which the C/G-type objects 
accreted from anhydrous minerals, organics and water ice. 
The recent discovery of comets orbiting among the main-belt 
asteroids, at semi-major axes similar to Ceres' [Hsieh & Jewitt, 
2006] support the possibility that small bodies can preserve part 
of their ices within the snow line region, which is defined as 
the heliocentric distance for which the temperature equals the 
condensation temperature of water (~5 AU). 

Other studies supporting a wet history for Ceres come from 
meteorites. Although no meteorites have been convincingly 
linked to Ceres [Feierberg et al., 1980; Jones et al., 1990; 
Sato et al., 1997], the C-type asteroids (whom G-type is a 
subclass) display a low albedo, and are thought to be the source 
of the Carbonaceous Chondrites (CC) meteorites. It has been 
shown that water ice could be stable inside CC meteorites 
over 4.5 Gyr [Fanale & Salvail, 1989, and references therein], 
thus supporting the possibility that the Main Belt has likely 
experienced a hydrated stage in its history (at least its outermost 
part). Moreover, spectral studies of Ceres reveal a strong 



3.07 jt/m absorption band characteristic of hydrated minerals 
[Jones et al., 1990; Sato et al., 1997]. This particular band has 
been the center of interest of many studies: In the early 80's, 
Feierberg et al. [1980] and Lebofsky et al. [1981] associated it 
with the signature of water frost on Ceres surface and predicted 
the possible existence of a polar cap. A decade later, King et al. 
[1992] fitted the 3 yum absorption with saponite, which is an 
ammonium-bearing phyllosilicate whose presence in CV and 
CI meteorites have been suggested by Zolensky & McSween Jr. 
[1988]. Later, Vernazza et al. [2005] found that crystalline water 
ice mixed with ion-irradiated asphaltite could reproduce a better 
fit to this feature. Recently, Rivkin et al. [2006] reviewed this 
past work and found that an hydrated iron-bearing phyllosilicate 
identified as cronstedtite, plus a few percent of carbonates 
could also fit adequately the near-infrared spectrum of Ceres. 
This recent interpretation could be supported by the detection 
in the mid-infrared range of emission features attributed to 
carbonates [from Cohen et al., 1998]. In addition, iron-bearing 
minerals have also been invoked by Feierberg et al. [1980] 
and Vilas & Gaffey [1989] to explain other spectral features 
present in the visible and near-infrared spectra of Ceres: 0.4 //m 
cutoff, 0.60 fim and 0.67 fim weak absorption bands and 1 yum 
shallow band [see also Vilas & McFadden, 1992]. Given the 
radical different interpretations of the 3.07 /mi band reported 
over the past decade, it is quite impossible to conclude on its 
exact nature. A persistent conclusion though, consists to support 
the presence of hydrated minerals, or residues from aqueous 
alteration, in the regolith material. 

As discussed above, surface water ice is not stable at 
distances smaller than 5 AU and is expected to sublimate if 
exposed directly to solar radiation [Fanale & Salvail, 1989]. 
Indeed, water ice migrating from the mantle region could pos- 
sibly reach the surface but will escape on very short time-scale 
[Nazzario et al., 2003]. As a result, direct detection of water 
ice vaporization from Ceres surface might be possible from 
the surroundings of a fresh impact crater, or from cracks in the 
sub-surface layers. Water escaping from Ceres polar region has 
possibly been detected in the early 90's [A'Hearn & Feldman, 
1992], although this 2-cr detection could never be unambigu- 
ously confirmed. Another supporting element comes from 
the relaxed shape of Ceres, which could be explained by the 
presence of large amount of ice in its interior [McCord & Sotin, 
2005]. The measurements of J2 obtained from HST and Keck 
both clearly support a model of a differentiated Ceres with a 
volatile-rich mantle, rather than a homogeneous interior model 
(see 3.3). 



4.2. Near-infrared maps 

To better represent the distribution and spatial extent of the 
features observed on the surface of Ceres (albedo, geological 
marks), we projected our high-angular resolution images into 
maps. The following sections are used to describe in details the 
various steps involved in the process of map projection, and the 
subsequent analysis of Ceres surface properties. 

Geometry: Any projection of an ellipsoidal shape onto a 
plane introduces deformations [Greeley & Batson, 1990]. We at- 
tempted to minimize these effects by choosing the cylindrical 
equidistant projection, which maps the surface of the asteroid 
onto a cylinder tangent to its equator and conserves the distances 
along the meridians. As a result, this projection minimizes the 
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Fig. 4. J-, H- and K- band maps of Ceres covering ~80% of Ceres' surface. The areas in white are terra incognita due to the limited 
ROI and rotational phases not imaged from Keck (see text). The theoretical resolution elements at J/H/K bands are shown at the 
bottom-right corners. Albedo variations are within ±6% around the mean surface value for each map. We estimate the errors to be 
2.5% maximum (see Fig. 5). The color scale is common for the three filters. Several round shaped features are visible including 
a dark feature with a bright center spot at ~(234°, +23°) (named "A") and a dark region at ~(125°, +20°) named "Piazzi" by 
Parker et al. [2002] (here "B"). Two other dark features are visible in the Southern hemisphere at ~(80°, -20°) and ~(285°, -35°) as 
well as other smaller features elsewhere (see Fig. 7). Any feature or albedo distribution present in the three maps has a very low 
probablity to be an artifact, with the exception of the diagonal stripes seen across the surface in H- and J-band (see text and Fig. 
5). One hemisphere (0° to 180°) appears to be sensibly brighter (1%) than the other at these wavelengths. A dark region running 
North/South is located at the boundary between the two hemispheres and is present in all maps. 



deformations of Ceres' equatorial area, which is seen at high- trie latitudes suffer stronger deformations after projection, but 
est resolution thanks to its small obliquity. Higher planetocen- 
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the impact is mitigated since these regions correspond to areas 
imaged tangentially, at a lower equivalent spatial resolution. 

Region of interest: We produced albedo maps of Ceres using 
the pixels located within an ellipsoidal Region Of Interest (ROI) 
centered on the image, and whose semi-axes were equal to 80% 
of the corresponding projected semi-axes on Ceres. The ROI cor- 
responds to 64% of the projected surface for each image (see 
Fig. 6). Ignoring the pixels near the edge of Ceres was based on 
several considerations: 

1- Although all images used to produce the maps were care- 
fully cleaned before deconvolution (including from corre- 
lated noise), some of them still presented artifacts after 
restoration of their optimal resolution, particularly near dis- 
continuities such as the limb and terminator. 

2- It appeared difficult in a few cases to restore optimally both 
the contour of the object and the surface details, even if 
Mistral is optimized to minimize the "ringing" effect intro- 
duced by the deconvolution of sharp edge objects. 

3- Finally, the resolution per pixel being highest at the center of 
the disk of Ceres, the use of the pixels located near the edge 
would strongly degrade the optimal resolution of our final 
product. 

Phase angle correction: As seen on Fig. 1, a phase angle of 
just a few degrees at near-infrared wavelengths produces a strong 
gradient of luminosity across the disk of Ceres. This effect had 
to be corrected prior to combine our individual images onto 
single maps. Several diffusion laws were investigated to model 
it, such as adopting a simple linear gradient, or using more 
complete models such as provided by the Lambert, Lommel- 
Seelinger, Minnaert and Hakpe laws (including single-scattering 
and multiple-scattering effects) [see Hapke, 1993]. The linear 
gradient was adopted since it provided the best fit to Ceres sur- 
face (residuals produced were nearly twice smaller than in the 
case of the Hapke model). 

Projection: The main difficulty to project an image of Ceres 
into its planetocentric referential (defined by its planetocen- 
tric longitude (A) and latitude (<p) as recommended by the IAU 
[Seidelmann et al., 2005]), resides in the accurate determination 
of its geometrical center. As the distribution of Ceres gravity 
field is unknown, we made the assumption that its center of 
mass, which is the center of the planetocentric referential, co- 
incides with its geometrical center. The conversion into plane- 
tocentric coordinates was based on our determination of Ceres 
dimensions (see 3.3) and the orientation of its spin axis (SVC, 
see 3.1), which is defined by the north pole angle (p„) and by the 
Sub-Earth Point coordinates (SEP^, SEP^). We used our value 
for Ceres rotation period (see 3.2) and the Eproc ephemeris gen- 
erator to obtain the SEPj, SEP^ and p„ at the time of each ob- 
servation. We then projected the images onto half-degree grid- 
ded maps to sample all areas of Ceres with sufficient resolution, 
from the equatorial region, up to the edge of the ROI. This fine 
grid map helped recover the smallest scale information from our 
set of 360 images. 

Combination of images into maps: No absolute photometric 
calibration was obtained for our data. We therefore used a near- 
infrared disk-integrated spectrum of Ceres [R.Binzel, personal 
communication] normalized to unity at 0.5 fim to calculate the 
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Filter 


Average Images 


Maximum Images 


Coverage 




(# per pixel) 


(# per pixel) 


(% total surface) 


J 


16 


63 


81.8 


H 


14 


54 


74.2 


K 


19 


72 


80.0 



Table 6. Average and maximum number of images used to pro- 
duce a single pixel of the final J/H/K bands maps, and their cor- 
responding surface coverage. 



equivalent disk-integrated photometric value for each band, and 
normalize our albedo maps with respect to each other. Prior 
to combine the different views into maps for each wavebands, 
we corrected the intensity variations of the individual projec- 
tions caused by the differential atmospheric absorption across 
our wavelength range, as well as the quality of the AO correc- 
tion. We first projected each image onto the planetocentric refer- 
ential of Ceres. Then we adjusted the brightness level of spatially 
adjacent maps by measuring the flux ratio over their overlap- 
ping area and applying the corresponding re-normalization coef- 
ficient. After correction, the maps could be combined to produce 
the albedo maps presented in Fig. 4 for each waveband. Each 
pixel of the final maps (covering a quarter of a square degree) 
was obtained by combining the corresponding pixels from the 
individual projections using a gaussian-weighted average func- 
tion (we chose a gaussian with a standard deviation of 5 pix- 
els). The largest weight (weight value = 1) was attributed to the 
pixel providing the best spatial resolution, while the pixels with 
a lower resolution were assigned a lower weight. Table 6 repre- 
sents the average and maximum number of images used to pro- 
duce a single pixel of the final maps for each waveband, as well 
as the percentage of surface coverage of Ceres (~80%). In or- 
der to facilitate comparison between our near-infrared data and 
the UV/Visible (223, 335 and 535 nm) HST maps, we chose to 
adopt the same reference meridian than published by [Li et al., 
2006a]. 

Maps description: The J-, H- and K-band maps shown in 
Fig. 4, and covering ~80% of Ceres' surface (see Table 6), are 
the result of combining 126, 99 and 135 individual projections 
respectively. We also derived error albedo maps (Fig. 5) by 
measuring, for each pixel, the intensity dispersion across the 
individual views. The theoretical size of the resolution element 
for J-, H- and K-band is 36.8 km, 47.4 km and 62.9 km respec- 
tively (corresponding to 4.4°, 5.6° and 7.5° at the equator). The 
major features sustain diameters of ~180 km (A and B) but 
smaller features can be seen in all three maps down to ~50 km 
scale. Although the theoretical resolution is highest in J-band, 
Fig. 4 shows that the final resolution is nearly equivalent across 
our three bands (i.e. ~60 km at equator). The degradation of 
the H- and J-band resolution is due to a more variable PSF at 
these wavelengths, which is supported by the larger photometric 
error (Fig. 5) derived for these bands. The amplitude of the 
albedo variation is within ±6% around the mean surface value 
for each band. The error maps show that the albedo maps (Fig. 
4) display an increasing error with decreasing wavelength: the 
1-cr uncertainty is smaller than 1% in the equatorial area in K 
band while it is estimated to be ~2% in J-band. The shape of 
the error distribution around 60° longitude reveal that the linear 
oblique features seen in J- and H-band are noise and should be 
ignored from our analysis. The K-band map displays the lowest 
noise level, due to the highest Strehl ratio delivered by the AO 
system at these wavelengths and should be considered as the 
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Fig. 6. Comparison between a synthetic view of Ceres (left) and 
a single deconvolved image (right). The overplotted dotted el- 
lipse corresponds to the ROI limit (80% of projected axis). The 
synthetic view appears "smoother" than the original view be- 
cause it is the result of the weighted average of several decon- 
volved images. The main albedo features are seen in both the 
synthetic and the original views, while other features did not re- 
sist the weighted average of several images, which acted as a 
filter to remove artifacts present in a single frame. 

any future ground-based, or space based observations of Ceres. 



Fig. 5. 1-cr deviation map expressed in percent with respect to 
the map mean value for each filter. The error is greater in J-band 
than in K-band as it could be expected from the lower strehl 
ratio at smaller wavelengths. The error also grows with latitude, 
which is due to the nearly equatorial view of Ceres as seen from 
an Earth observer. 



most accurate of the three. 

Some of the main albedo features in our near-infared maps 
can also be seen in the UV/Visible HST maps published by 
Li et al. [2006a], like the large bright area around (125°, +20°) 
and the dark spot at (130°, -24°). The fact that these features 
are visible in all wavelengths suggests that they are geological 
features like basins or impact craters... But the UV/Visible 
and near-infrared maps show also discrepancies: the bright 
feature at (115°, -30°) visible in the UV/Visible maps is not 
present in our maps, nor the dark feature located around (45°, 
+ 10°). Whereas all the discrepancies cannot be linked to surface 
properties (the dark annulus of the "A" feature (234°, +23°) 
is not visible in the UV/Visible map due to its size which is 
comparable to the resolution element provided by the HST), 
such variations in their spectral behavior suggest compositional 
differences between these regions. The next section will provide 
a discussion on possible composition and geological origin. 



4.3. 3-D model 

We re-projected our multicolor maps of Ceres onto its 3D shape 
model. Such model is used to check the distribution of albedo 
features present in our final maps against the features seen in 
our original images after deconvolution (see Fig. 6). This model 
can also be used to predict 1 what Ceres surface would look like 
at any epochs, which will certainly be useful when preparing 



1 http://www.imcce.fr/page.php ?nav=en/ephemerides/formulaire/form. 



4.4. Three-band spectra 

We investigate the nature of the major features seen in our maps 
by reporting their relative brightness variation at J-, H- and K- 
bands. We selected 10 type of area (6 bright and 4 dark, noted 
b, and d, respectively) and represented in Fig. 7 their photomet- 
ric measurements normalized to Ceres integrated spectrum. The 
photometric value obtained for each band is the result of aver- 
aging the photometric measurements over an area equivalent to 
a resolution element. We estimate the error on the relative pho- 
tometry as J^j ct 2 , where cr, is the standard deviation for each 

pixel as read in the error map (Fig. 5). The results are shown in 
Fig. 7 with their 3-cr error bars. 

Whereas spectral variation with the rotational phase of Ceres 
has never been reported, our data show differences of spectral 
behavior across the surface. The analysis of Fig. 7 reveals a 
clear trend of the bright features (left) to display a higher H-band 
albedo (with respect to J- and K-band) than the rest of Ceres sur- 
face. On the contrary, the dark features (right) do not display a 
similar trend. This may indicate a common origin for the bright 
features, whereas the dark regions may be the result of various 
surface processes, or represent different level of surface aging. 

We then compare the 3-band spectra to those of various labo- 
ratory compounds. We used the ASTER 2 and RELAB 3 spectral 
libraries, to obtain the equivalent near-infrared spectra of vari- 
ous compounds predicted to be present on Ceres. We report in 
Fig. 8 their broad band photometric values, normalized to Ceres 
disk-integrated spectrum, each mineral being mixed with a dark 
neutral compound, using the mixing ratios given in Table 7. 
Due to the low spectral resolution provided by our broad band 
imagery, it is not possible to identify unequivocally the com- 
pounds present on the surface of Ceres. We thus based our study 

2 http://speclib.jpl.nasa.gov 

3 http://www.planetary.brown.edu/relab/ 
ephephys.php 
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Fig. 7. Three-band spectra of selected surface features normal- 
ized to Ceres disk-integrated spectrum. The left (respectively 
right) panel shows the photometric points obtained for the bright 
(respectively dark) features. All spectra are shifted vertically by 
0.045 (respectively 0.065) to improve clarity. The letters posi- 
tioned left of each spectrum are used to locate the corresponding 
features on the map. The K-band resolution element is shown at 
the bottom-right corner. 



to the comparison between the possible surface compounds pro- 
posed by Rivkin et al. [2006], who showed that carbonates like 
Siderite, Dolomite or Calcite mixed with phyllosilicates pro- 
vide an excellent fit to the 3 /mi region of Ceres. Carbonates 
could also explain the mid-infrared spectral emission features 
detected by Cohen et al. [1998]. We also included an orthopy- 
roxene (Enstatite), a clinopyroxene (Augite) and Olivine, the 
most abundant elements in the Solar System. A water ice frost 
spectrum is also reported. 

Although the compositional information returned at such a 
low spectral resolution should be considered with caution, we 
can use the comparison between the spectral behavior measured 
on the surface and that of the laboratory samples to constrain 
and discuss the possible surface composition of Ceres. There are 
similarities between the spectral behavior of the bright features 
seen in Fig. 7 (left) and that of phyllosilicates and carbonates 
in Fig. 8. The Calcite and Montmorillonite "spectra" display the 
same shape as the majority of the bright features. For instance, 
the Montmorillonite reproduces quite well the behavior of the 
b4 feature as well as the bright region b2 surrounding d3. On 
the other hand, the spectral behavior of the Siderite (a iron-rich 
Calcite), Cronstedtite (an iron-bearing phyllosilicate), or Augite 
(a clinopyroxene) only match the spectrum of a small percent- 
age of Ceres surface, and for small mixing ratios show that all of 
them might be minor compounds of the regions discussed in this 
section. Although, igneous rocks like pyroxenes and olivine are 
not expected to be present on the surface of primitive asteroids, 
these compounds are likely minor components of the regions dis- 
cussed in this section if present at all. 




Montmorillonite 



Phyllosilicates 



Augite 




I ... I ... I ... I 



5? 



1.0 1.2 1.4 1.6 1.8 2.0 2.2 
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Fig. 8. Three-band spectra of selected carbonates [CO? - ], 
phyllosilicates [S1O3], pyroxenes [(Si,Al)20e], olivine 
[(Mg,Fe)2Si04], and water ice frost [H2O] normalized to 
Ceres disk-integrated spectrum. Error bars are negligible. 



Component 


Grain Size 


Mixing ratio 


Sample 




(jum) 


(%) 


No. 



Calcite" 


0-45 


100 


C-3A 


Dolomite " 


0-45 


50 


C-5A 


Siderite" 


0-45 


4 


C-9A 


Crondstedtite 


0-45 


5 


CR-EAC-021 


Montmorillonite " 


0-45 


50 


PS-2D 


Enstatite " 


0-45 


30 


IN-10B 


Augite " 


0-45 


2 


IN- 15 A 


Olivine ' 


0-50 


5 


DD-MDD-013 


Water Frost " 


-10 


2 


FROST.SNW 



Table 7. Grain size and mixing ratio for the selected compounds. 
The last column indicates the sample reference in the spectral 
libraries (a: ASTER and r. RELAB). 



The interpretation of the dark features is more complex. 
Indeed, the lower H-band value that characterizes half of the 
dark features (di, d2) does not match the behavior seen in phyl- 
losilicates and carbonates, which are the major compounds pre- 
dicted to be present over Ceres (see 4.1). Such a drop in H-band 
matches better the behavior of Enstatite, but its high density of 
~3 200 kg.irT 3 and the mid-infrared mismatch between Ceres 
and Enstatite spectra makes its presence unprobable. Another 
possible matching compound is water ice (Fig. 8). Water ice is 
not expected to be found in a stable form over Ceres, but "dirty" 
ice, i.e. ice or frost mixed with other minerals, could be present 
and reproduce the more elevated K-band point seen in the "spec- 
tra" of these features. Another point to consider is the relatively 
high planetocentric latitude of these dark features. None of them 
is found in the equatorial region, most being located above 30- 
40° latitude. If dirty ice exists on the surface or sub-surface lay- 
ers of Ceres [McCord & Sotin, 2005; Mousis & Alibert, 2005], 
it would be expected to be more stable at higher latitudes, where 
surface temperature is lower. Nevertheless, we cannot ascertain 
from these data alone that ice is present on Ceres. Higher spectral 
resolution, coupled with the high-contrast and spatial resolution 
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provided by adaptive optics, is required to investigate in detail 
the composition of the main features seen on its surface. The 
limits of such a broad-band analysis in constraining the compo- 
sition of the main features is illustrated by the bright feature bi, 
which appears to match quite well the behavior expected for wa- 
ter frost, but it is located at low latitude and differ in reflectivity 
from the dark features discussed above. 

The two main observed features A (b4, d4) and B (d3 and 
the surrounding region b2), were both referenced as #5 and #2 
by Li et al. [2006a]. These two features are remarkable because 
whereas they sustain large physical dimensions (©a ~ 180 km 
and ©g ~ 350 km), their spectral behavior appears to be homo- 
geneous over such large area, which might point to a same com- 
position and/or resurfacing history. The external annulus of the 
B feature (b2) shows the same uniformity. This annulus shape is 
reminiscent of a large cratering event. The bright central regions 
of feature A resembles the central peak seen in craters originat- 
ing from high energetic impact, which would then be subject to 
different aging processes than the lower altitude neighboring ar- 
eas. At smaller physical scale, similar differences of albedo have 
been reported on asteroid 25143 Itokawa between the central 
part and surrounding areas of small craters, this time via grain 
sorting [Miyamoto et al., 2007]. 

In summary, our AO study permits to map the albedo varia- 
tions over the surface of Ceres down to a 40-60 km scale and 
investigate whether these variations correspond to changes in 
the composition. Also, the high-spatial resolution capabilities re- 
turned by AO make possible to search for areas of distinctive 
signature (e.g. icy rich spots), which would remain undetected 
otherwise in disk-averaged studies of Ceres. A dusty regolith a 
few centimeters thick [Webster et al., 1988], created by the im- 
pact of micrometeorites and possibly larger bodies, has been pro- 
posed [Witteborn et al., 1999; Lim et al., 2005] as a possible ex- 
planation for the shallow spectral signatures in the visible and 
near-infrared spectrum of Ceres and its small albedo constrast. 
But some small areas visible in our high-angular resolution im- 
ages of Ceres could correspond to places on the surface where 
the regolith material has been cleared by "recent" endogenic or 
exogenic activity to expose more pristine material from the sub- 
surface layers. 

5. Conclusions 

We present the results of high-angular resolution near-infrared 
(~[1. 1,2.3] jum) imaging observations conducted with the adap- 
tive optics system at the Keck II telescope in Hawaii. Imaging 
Ceres over nearly one rotation period allows us to derive its main 
physical characteristics: rotation period (P s = 9.074 10^J5oooi4 
h), spin vector orientation (EQJ2000.0 a = 288° ± 5°, '8 = 
+66° + 5°), size ( a = b = 479.7 + 2.3 km, c = 444.4 ± 2.1 km ) 
and density (p - 2 206 ± 43 kg.trT 3 ) and produce the first near- 
infrared maps of its surface, revealing a wealth of surface fea- 
tures (also seen in the UV/Visible wavelengths [Li et al., 2006a]. 
We discuss the nature and composition of the main albedo marks 
present on its surface. We do not detect topographic features over 
18 km high, which agrees with predictions for a body of this size 
whose internal composition could still contain a large amount 
of water ice. By comparison with laboratory samples, we found 
Ceres surface composition to be consistent with phyllosilicates 
such as Montmorillonite, and carbonates such as Dolomite or 
Calcite. Some dark regions above 30-40° latitude display the 
characteristics expected from water-frost mixed with minerals. 
Unfortunately, broad band imaging does not provide sufficient 
spectral resolution to investigate further the composition of dis- 



tinct geological features present on Ceres, and future spectro- 
scopic observations combined with high-angular resolution are 
required to better constrain its mineralogy. 
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